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ADJOINT-BASED GRADIENT DRIVEN METHOD FOR IDENTIFYING 
UNKNOWN PARAMETERS OF NON-LINEAR SYSTEM MODELS 

CROSS-REFERENCES TO RELATED APPLICATIONS 

[0001] This application claims the benefit of U.S. Provisional Application 
No. 60/455,083, filed March 13, 2003. 

FIELD OF THE INVENTION 
[0002] The present invention generally relates to a method for dynamic 
system parameter identification and, more particularly, to an adjoint-based 
gradient driven method for identifying nonlinear system dynamic parameters 
for automotive powertrain systems and subsystems. 

BACKGROUND OF THE INVENTION 
[0003] Controllers for systems and subsystems, such as engine powertrain 
systems and subsystems, are increasingly being implemented using modern 
control techniques. Modern control techniques are derived from mathematical 
models that mathematically describe the dynamic behavior of the system to be 
controlled. The mathematical model of such systems can be derived, for 
example, from one or more physical laws, or from experimental data using 
known regression techniques. 

[0004] No matter how the mathematical system model is derived, it should 
be verified to ensure that it represents the real system with sufficient accuracy. 
To do so, the computed outputs of the system model may be compared to 
actual, experimentally sensed data that describe the same output. 
Mathematical methodologies, such as least square fit algorithms, may then be 
used to provide some measure of difference between the computed and actual 
outputs. If this indicates that the fit is not acceptable, appropriate 
modifications may be made to the model, to the initial conditions, and/or to 
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other characteristics, to improve the fit. This process is iterative in nature, and 
is conducted until the difference between the computed and actual outputs 
sufficiently converge, which can be both time consuming and costly. 
[0005] Accordingly, it is desirable to provide a method for identifying the 
5 unknown parameters of a non-linear dynamic system that is computationally 
efficient as compared to presently known methods, is relatively less time 
consuming than presently known methods and, thus, is less costly than 
presently known methods. Furthermore, other desirable features and 
characteristics of the present invention will become apparent from the 
10 subsequent detailed description of the invention and the appended claims, 
taken in conjunction with the accompanying drawings and this background of 
the invention. 



BRIEF SUMMARY OF THE INVENTION 
15 [0006] A method is provided for identifying the unknown parameters of a 
non-linear dynamic system model that has one or more system inputs. The 
method includes determining a governing state equation from the system 
model. A generalized cost function that represents a performance objective 
for the system is determined. An adjoint equation is determined based at least 
20 in part on the governing state equation. A gradient is determined based at 
least in part on the adjoint equation. The governing state equation, the adjoint 
equation, and the generalized cost function are supplied to a processor. The 
processor is then caused to iteratively determine changes in the generalized 
cost function that result from incremental changes in arbitrarily chosen values 
25 of one or more of the unknown model parameters to thereby identify the 
unknown model parameters. 

[0007] A computer-readable medium containing computer executable code 
that instructs a computer to perform the above-described method is also 
provided. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0008] The present invention will hereinafter be described in conjunction 
with the following drawing figures, wherein like numerals denote like 
elements, and: 

5 [0009] FIG. 1 is a flowchart depicting a computational method for 
identifying unknown model parameters of a non-linear dynamic system model, 
such as a model of the torque converter of FIG. 2, according to an embodiment 
of the present invention; 

[0010] FIG. 2 is a simplified functional block diagram of an automobile 
10 powertrain system; 

[0011] FIG. 3 is a simplified cross-section of a torque converter that may be 
included in the powertrain system of FIG. 2; 

[0012] FIG. 4 is a simplified representation of a control volume based on the 
torque converter depicted in FIG. 3; and 
15 [0013] FIG. 5 is a graph depicting the response of a non-linear system model 
of the torque converter of FIGS. 3 and 4, in which the model parameters were 
determined using the process illustrated in FIG. 1. 



20 DETAILED DESCRIPTION OF THE INVENTION 

[0014] The following detailed description of the invention is merely 
exemplary in nature and is not intended to limit the invention or the 
application and uses of the invention. Furthermore, there is no intention to be 
bound by any theory presented in the preceding background of the invention 

25 or the following detailed description of the drawings. 

[0015] A flowchart depicting a generalized process 100 for determining the 
system parameters of a non-linear dynamic system model according to an 
exemplary embodiment is shown in FIG. 1. It should be appreciated that the 
process 100 is preferably implemented, in whole or in part, on an 

30 appropriately programmed general purpose computer. Alternatively, a 
specialized device may be designed and constructed to implement the method 
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in hardware, firmware, software, or combination thereof. In the following 
description, the parenthetical references in FIG. 1 correspond to the particular 
reference numerals of the methodological flow illustrated therein. 
[0016] As FIG. 1 indicates, the first step to identify system model 
5 parameters is to determine a governing state equation from the system model. 
An general form for the governing state equation for the nonlinear systems we 
may consider with this approach is: 

q = N(q,9,\|/) on 0 < t < T, for q = q o , at / = 0, (1) 

10 

where 9 is the parameter vector being sought, and N is a continuous, 
differentiable nonlinear function of the state vector q, the parameter vector 0 , 
and the exogenous input vector \|/. 

[0017] After the state equation is determined, an adjoint optimization 
15 procedure is defining is then implemented to identify the model parameters 
(104-110). Adjoint analysis is a mathematical tool that is used to determine 
the gradient information that is central to efficient high-dimensional 
optimization strategies. Although adjoint analysis has been used previously to 
conduct control input optimization, it has not been used to conduct model 
20 parameter identification, which is the subject of the present invention. As 
FIG. 1 indicates, the first step in the adjoint optimization procedure is defining 
a cost function that represents the performance objective for the process (104). 
A generalized expression of the cost function (J) is shown below. 

(2) 

In the context of the present embodiment, the objective is output validation of 
the model to ensure the parameters were correctly determined. The error 
30 between the model output Cq and the measured data y is penalized by 
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adjusting the weighting matrix Q. Moreover, since the measured data 
includes noise, the single sample for the initial condition q o cannot be 
provided with 100% certainty. Hence, the model initial conditions are 
permitted to vary by adjusting the weighting matrix R 0 . The parameter values 

5 can be constrained near reference values 8 by adjusting the cost weighting 
matrix R . 

[0018] The model output validation objective is met when equation (2) is 
minimized with respect to the parameters 0 being sought, subject to the 
governing state equation (1). Knowing the effect of parameter changes on 
10 equation (2) helps search for this minimum. Now, letting q' be the small 
change in the state q when the parameters are perturbed a small amount, and 
defining the following equation: 

Lq' A(d/3r-A)q\ (3) 

15 

which is the linearization about a trajectory q(0 ,q o ), a perturbation equation 
can be written as follows: 

L q' = B e Q on 0 < t < T, for q' = 0, at t - 0, (4) 

20 

where 0 is a parameter perturbation from 0 , that drives the state perturbation, 
q\ which has an associated non-negligible perturbation cost function defined 

as: 

25 J = {"(Cq-y) QCq'dt + (C 0 q 0 - yJ*tf 0 C 0 q' o +{b-o)*R»\ (5) 

[0019] The sensitivity of equation (5), the perturbation cost function, with 
respect to the parameter perturbations 0 may then be found, since this 
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determines how to minimize the cost function. To do so, the adjoint analysis 
proceeds by introducing an inner product defined as follows: 

<r, q' > = £r*q'rf/, (6) 

5 where the asterisk implies the transpose operation. This inner product yields 
the adjoint identity: 

<r, L q* > = <L*r, q' > + b, (7) 

10 where r is introduced as the adjoint state, such that L*r = {-dldi - A*)r, and 
such that b = [r*q f ] t =T - [r*q'] t =o. From the adjoint identity, and choosing the 
adjoint state equation to be driven by the model validation error defined in the 
cost function, the following adjoint equation is obtained: 

15 £*r= C*Q(Cq-y) on 0 < t < T, for r = 0, at t = T. (8) 

[0020] Once the adjoint equation is determined (106), a perturbation cost 
function is then easily computable. Specifically, by combining equations (4) 
(the perturbed state equation), (7) (the adjoint identity) and equation (8) (the 
20 adjoint equation), the following alternative expression for the perturbation cost 
function (5) may be obtained: 

3 = [fo B ° Tdt + *(° " 5 )] °' + fo*. - y o ) + r(0)]V„ (9) 

25 [0021] The unknown model parameters are then identified using equation 
(9) by iteratively determining changes to the cost function (equation (2)) with 
respect to changes in the system initial conditions and the model parameters. 
Specifically, arbitrary values for both the initial conditions q o and the 
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unknown model parameters 0' are selected, and the cost sensitivities 
associated with these changes are determined. In particular, the following 
gradients are determined (108), and then iterative calculations are conducted 
(1 10) using these gradients to determine the sensitivities and the minimization 



[0022] A particular preferred embodiment of the iterative calculations that 
10 are carried out as part of the generalized functional step 1 10 are also illustrated 
in FIG. 1, as steps 111-125. The skilled artisan will readily appreciate and 
understand the steps illustrated in this portion of the flowchart, and will 
therefore not be described herein in further detail. 

[0023] The determined model parameters are then used to realize the 
15 original system model state equation using the same state space notation that is 
used in the adjoint analysis process (112). Thereafter, the realized system 
model state equation may be validated against one or more sets of 
experimental data (114). 

[0024] An example of how the above process can be applied to an 
20 automobile powertrain system will now be discussed. Before doing so, an 
exemplary powertrain system will first be described. Turning first to FIG. 2, a 
schematic diagram of an automobile powertrain system 200 is depicted. The 
powertrain system 200 includes an engine 202 and a transmission 204. The 
engine 202 is the prime mover of the vehicle into which the powertrain system 
25 200 is installed. The engine 202 responsive to driver input from a throttle 
pedal 203 to a powertrain controller 205, and generates the torque necessary to 
accelerate the vehicle to a desired velocity, and to maintain the vehicle at this 
desired velocity. The torque generated by the engine 202 is supplied, via an 
engine flywheel 206, to the transmission 204. The transmission 204 in turn 
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= CX(a-yJ + r(0). 
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couples the torque supplied from the engine 202 to various numbers of driven 
wheels 207 via selected ones of a plurality of fixed gear ratios, which are 
housed within a transmission gearbox 208. The transmission 204 additionally 
includes a torque converter 210, which provides a hydrodynamic coupling 
5 between the engine 202 and the transmission 204. A simplified cross-section 
of an exemplary embodiment of the torque converter 210 is illustrated in FIG. 
3, and will now be described. 

[0025] The torque converter 210, which is exemplary of a general torque 
converter that may be used in any one of numerous powertrain systems 200, 

10 includes a housing 302, a pump 304, and a turbine 306. The housing 302 is 
coupled to the engine flywheel 206. Thus, as the engine 202 rotates the engine 
flywheel 206, the housing 302 is also rotated at the same rotational speed. The 
pump 304 is a centrifugal -type pump having an impeller with a plurality of 
fins 308. The fins 308 are coupled to the housing 302, and therefore rotate at 

15 the same rotational speed as the engine 202. As the pump 304 is rotated, 
hydraulic fluid within the housing 302, which is preferably automatic 
transmission fluid, is thrown outwardly by the impeller fins 308 toward the 
housing 302. This creates a vacuum in the center of the impeller that draws 
more fluid into the pump 304. 

20 [0026] The turbine 306 includes a plurality of blades 310. The fluid that 
exits the pump 304 strikes the blades 310, causing the turbine 302, and thus 
the transmission 104, to rotate. The blades 310 are preferably curved, so that 
the fluid that enters the turbine 306 changes direction before it exits the center 
of the turbine 306. This directional change is what causes the turbine 306 to 

25 spin. As the fluid exits the center of the turbine 306, it is moving in a different 
direction than when it entered, which is a direction that is opposite that which 
the pump 304 is turning. Thus, the torque converter 210 may additionally 
include a stator 312. The stator 312 is centrally disposed between the pump 
304 and the turbine 306, and redirects the fluid returning from the turbine 306 

30 before it reaches the pump 302. 
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[0027] To apply the process 100 described above to the powertrain system 
200, and more particularly to the torque converter 210 subsystem, a model of 
the torque converter is first developed. To do so, the torque converter 210 is 
modeled, as shown in FIG. 4, as a combination of two fixed control volumes 
5 402 and 404, with the stator volume excluded. The first control volume 402 
contains the automatic transmission fluid in the pump 304, and the second 
control volume 404 contains the automatic transmission fluid in the turbine 
306. The first 402 and second 404 control volumes share a boundary 406 at 
the interface where the fluid leaves the pump 304 and enters the turbine 306, 
10 and vice-versa. The conservation of momentum equation for a fixed control 
volume, delineated below, is then applied twice to obtain two equations, one 
for each control volume: 

rxF s +l cv rxgpdV+T shafi = j-j^rxVpdV + ^rxVpV - dA 9 (12) 

15 

where V is the fluid velocity, p is the fluid density, CS is the control surface 
around each control volume 402, 404, r is the appropriate moment arm; F s is 

any force affecting the control surface; JV x gpdV is any force caused by 

gravity (g) that affects the control volume; and T^ft is the externally applied 

20 mechanical torque. 

[0028] Due to the symmetry of the torque converter 210, the net 
gravitational force (g) for both the pump 304 and the turbine 302 is zero. 
There are, however, three forces that affect the control surfaces. The first is 
skin friction loss along the path of the pump impeller fins 308 and the turbine 

25 blades 310. As is generally known, skin friction loss is a linear function of the 
fluid speed. The second force is the shear loss that is incurred at the boundary 
406 between the first 402 and second 404 control volumes. The shear loss is 
a linear function of slip (AN=N engine - N turbine ). And, the third force that acts on 
the control surfaces is the head loss that is also incurred at the boundary 406 

30 between the first 402 and second 404 control volumes. The head loss is a 
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generally known quadratic function of slip. As regards the externally applied 
mechanical torque, for the pump 304 this is any brake torque that is applied, 
and for the turbine 306 this the vehicle load torque transferred by the gear 
shifting device in the gear box 208. 
5 [0029] Now, by combining equation (12) for both of the control volumes 
402 and 404, and consolidating various terms in a generally well-known 
fashion, the following coupled expression for the torque converter system is 
obtained: 



•j^q @aY turbine + ^b^ turbine + turbine^ engine + engine ~ ^e^ engine + @f ^ turbine + ^g^ engine + turbine 
engine + turbine + turbine ^ engine + engine ~ engine + ^n^ turbine + engine + ^p^ turbine 

(13) 



[0030] For the purposes of system identification, the system output is chosen 
to be the full state, consisting of both engine speed ( N engine ) and turbine speed 

15 turbine)* anc * f° r simplicity, consider the engine torque (Y engine ), the turbine 

torque ( y turbine ), and the transmission shaft speed ( N output ), to be exogenous 

inputs, that is given signals or measurements for all time. However, the lack 
of available turbine torque measurements required some analysis, which show 
that the turbine torque can instead be approximated by: 
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N . 

v ~6 N +0v +0 engme Y (14) 

/turbine q output w r f output 1 w j ^ 'engine* V A ^V 

turbine 



where transmission shaft torque ( y oulput ) and transmission shaft speed ( N ouJpul ), 

are accessible measurements. Thus, rather than turbine torque, the 
25 transmission shaft torque is considered an exogenous input along with the 
engine torque and transmission shaft speed ( N omput ),. 

[0031] Now, rearranging equation (13) so that the time derivatives are on 
the left, appropriately substituting in equation (14), and consolidating constant 
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parameters, the governing state equation for the torque converter 210 is 
obtained: 



turbine 



4* *14 



&} N turbine + &4 N engine 
[&i2 N turbine + &[* N engine 



e 5 N engine turbine 
9\ 4 N engine J J |_ ^engine 
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Yengine 



(15) 



[0032] Turning now to the remainder of the process 100, the torque 
converter governing state equation (equation (15)) is applied with the adjoint 
equation (equation (8)), and the cost function (equation (9)). Then, iteratively 
10 calculating the gradients(equations (10) and (11)), the realization of equation 
(15), using the system model notation delineated in state equation (1) above is: 



N 



turbine 



N 



V = 



engine 



Y engine 
y ' output 

N 

output 



(16) 



15 It is noted that the over bars indicated in equations (16) above indicate that the 
variables are normalized. Specifically, the speed variables are normalized to 
the maximum engine speed, and the torque variables are normalized to 
maximum transmission shaft torque. 

[0033] The state space model (equation (15)) with the identified system 
20 parameters (dj.is) was validated against experimental data associated with a 1- 
2 upshift during a 30% constant throttle pedal maneuver. For this single set of 
data, no reference parameters are known, and therefore the parameter weight 
(R ) is zero. In addition, the weight on the output error (Q) is chosen as the 
identity matrix, and a non-zero value of the initial condition weight (R 0 ) is 

25 chosen. The results are shown in FIG. 5, which depicts normalized turbine 
speed and normalized engine speed versus time for both the simulation model 
output 502 and 504, respectively, and the experimental data 506 and 508, 
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respectively. Values of measured engine torque, transmission shaft torque, 
and transmission shaft speed were used as system inputs. 

[0034] It will be appreciated that the identified model should be validated 
against more than one data set, since experimental measurements may vary for 
5 a specific operational maneuver, such as the 1-2 upshift during a 30% throttle 
pedal maneuver described above, from experiment to experiment. Thus, to 
provide increased assurance of model validity, the parameter identification 
process 100 should be validated against an ensemble of data sets. For 
example, in a particular embodiment, the torque converter system model 

10 parameters (Oj-js) were validated against an ensemble of 10 data sets. In 
conducting such an extended validation, the parameters and initial conditions 
found with the single data set were used as the initial guesses. The cost 
function now includes the integral of the error from all 10 data sets over the 
time interval [0, 7], as well as the variation in initial conditions with respect to 

15 the 10 different measured initial conditions. 

[0035] When a plurality of data sets is used, the overall process 100 is 
unchanged from that used for a single data set; however, the process 100 now 
has information available from all of the data sets in the ensemble. Therefore, 
the model is more generally validated. Indeed, it was found that the torque 

20 converter model that was validated based on the ensemble of data sets 
provided a response that closely approximated a 1-2 upshift for an arbitrary 
constant pedal maneuver. 

[0036] The above-described methodology provides for the identification of 
the unknown parameters of a non-linear dynamic system in a manner that is 

25 computationally efficient as compared to presently known methods. The 
method is relatively less time consuming than presently known methods and, 
thus, is less costly than presently known methods. Moreover, the method 
provides the flexibility to allow the initial conditions used during the 
parameter identification process to be varied from preferred values, the 

30 method allows the state values to stay closer to nominal values, which ensures 
physically meaningful results are provided, and allows more flexibility. 
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[0037] While an exemplary embodiment(s) has been presented in the 
foregoing detailed description of the invention, it should be appreciated that a 
vast number of variations exist. It should also be appreciated that these 
exemplary embodiments are only examples, and are not intended to limit the 
scope, applicability, or configuration of the invention in any way. Rather, the 
foregoing detailed description will provide those skilled in the art with a 
convenient road map for implementing a preferred embodiment of the 
invention. It being understood that various changes may be made in the 
function and arrangement of elements described in an exemplary preferred 
embodiment without departing from the spirit and scope of the invention as set 
forth in the appended claims. 



